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Abstract. The classical Heun equation has the form 

Q{Z) J? + P(z) fz + V(Z) } S(Z) = °' 
where Q(z) is a cubic complex polynomial, P(z) is a polynomial of degree at 
most 2 and V(z) is at most linear. In the second half of the nineteenth century 
E. Heine and T. Stieltjes in [5], 1131 initiated the study of the set of all V(z) 
for which the above equation has a polynomial solution S(z) of a given degree 
n. The main goal of the present paper is to study the union of the roots of 
the latter set of V(z)'s when n — » oo. We formulate an intriguing conjecture 
of K. Takemura describing the limiting set and give a substantial amount of 
additional information obtained using some technique developed in [7], 



1. Introduction and Main Results 
A generalized Lame equation is a second order differential equation of the form 

d 2 d 



Q(*)^ + r(z) Tz + v( z )^( z ) = o, (i) 

where Q(z) is a complex polynomial of degree I and P(z) is a complex polynomial of 
degree at most l — l, see [T7]- It was first shown by Heine [S] that if the coefficients 
of Q(z) and P(z) are algebraically independent, i.e. do not satisfy any algebraic 
equation with integer coefficients then for an arbitrary positive integer n there are 
exactly (" + ~ 2 ) polynomials V{z) such that |T]) has a solution S(z) which is a 
polynomial of degree n. As was recently shown in [Uj for any equation |T]) with 
degQ(z) = l,dcgP(z) < I — 1, and any positive n the set Q3„ of all V(z) giving a 
polynomial solution S(z) of degree n is always finite and its cardinality is at most 
( n+l - 2 ). Below we concentrate on the classical case I = degQ(z) = 3 which is 
better known under the name the Heun differential equation, see e.g. [5] and study 
the union of all roots of polynomials V(z) belonging to 5J„ as n — > oo. Note that 
if I = degQ(z) = 3 then V(z) is at most linear and that for a given value of the 
positive integer n there are at most n+l such polynomials. 

No essential results in this direction seems to be known. One of the few excep- 
tions is a classical proposition of Polya, |10j claiming that if the rational function 
has all positive residues then any root of any V(z) as above and of any S(z) 
as above lie within Couvq where Couvq is the convex hull of the set of all roots of 
Q(z). 

Before we move further let us formulate appropriate versions of two main results 
of [TT] generalizing the above statements of Heine and Polya. 

Theorem 1. For any polynomial Q(z) of degree I and any polynomial P(z) of 
degree at most I — 1 
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• there exists N such that for any n > N there exist exactly (™ + ^ -2 ) polyno- 
mials V(z), deg V{z) = 1 — 2 counted with appropriate multiplicity such that 
|T]) has a polynomial solution S(z) of degree exactly n; 

• for any e > there exists N e such that for any n > N e any root of any 
above V(z) and S(z) lie in the e-neighborhood of Couvq. 

Applying the latter result to the situation I = 3, i.e to the Heun equation we can 
introduce the set V n consisting of polynomials V(z) giving a polynomial solution 
S(z) of ([1]) of degree n; each such V(z) appearing the number of times equal to its 
multiplicity. Then by the above results the set V n will contain exactly n + 1 linear 
polynomials for all sufficiently large n. It will be convenient to introduce a sequence 
{Sp n (X)} of spectral polynomials where the n-th spectral polynomial is defined by 

71+1 

Sp n (X) = JJ(A - t n! j), 
3=1 

where t n j is the unique root of the j-th polynomial in V„ in any fixed ordering. 
(Spn(X) will be well-defined for all sufheicntly large n.) 
Associate to Sp n (X) the finite measure 

^ n+1 

3=1 

where 5(z — a) is the Dirac measure supported at a. The measure fx n obtained in 
this way is clearly a real probability measure which one usually refers to as the 
root- counting measure of the polynomial Sp n (X). 

The starting point of this project was some numerical results for the distribution 
of roots of Sp n (X) obtained by the first author about 5 years ago and illustrated 
on the next figure. 
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Figure 1. The roots of the spectral polynomial Sp5o(X) for the 
classical Lame equation \Q{z)-^ + \Q'{z)-^ + V(z) j S(z) = 0, 
with Q(z) = z(z - 1) (z + | - i). 

Extensive numerical experiments strongly suggest that the following holds. 

Conjecture 1 (Shapiro- Tater). For any equation (TT]) the sequence {fi n } of the root- 
counting measures of its spectral polynomials converges to a probability measure \x 
supported on the union of three curved segments located inside Couvq and connect- 
ing the three roots of Q(z) with a certain interior point, see Fig. [H Moreover, the 
limiting measure [i depends only on Q(z), i.e. is independent of P{z). 

An elegant description of the support of \x was suggested to us by Professor 
K. Takemura, [T5] , 
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Denote the three roots of Q{z) by oi, 02, 0,3. For i <E {1, 2, 3} consider the curve 
7i given as the set of all b satisfying the relation: 

f ak l~ b-t 

/ \ 1 77 77 rdteR, (2) 

here j and are the remaining two indices in {1,2,3} in any order and the inte- 
gration is taken over the straight interval connecting a,j and Ofc. One can see that 
cii belong to 7^ and that these three curves connect the corresponding a, with a 
common point within Couvq. Take a segment of 7^ connecting Oj with the common 
intersection point of all 7's. Let us denote the union of these three segments by 

Conjecture 2 (Takemura). TVie support of the limiting root-counting measure fi 
coincides with the above Tq. 

The above description of Tq led us to the following reformulation of Takcmura's 
conjecture. 

Proposition 1. The above set Tq coincides with the continuum of minimal loga- 
rithmic capacity connecting the roots ofQ(z). 

Notice that Goluzin's classical problem of finding the continuum of minimal 
capacity connecting a given n-tuple of points in C was completely solved for n = 3 
by G. Kuzmina in [5], see also [5]. 

In the joint with Professor Takemura follow-up of the present paper |12j we will 
completely settle the above Conjecture [5] and Proposition [T] using some methods 
and results presented below. In the present paper generalizing the technique of [7] 
we study a different probability measure which is easily described and from which 
the measure fi (if it exists) is obtained by the inverse balayage, i.e. the support of fi 
will be contained in the support of the measure which we construct and they have 
the same logarithmic potential outside the support of the latter one. This measure 
will be uniquely determined by the choice of a root of Q(z) and thus we are in fact 
constructing three different measures having the same measure fi as their inverse 
balayage. 

1.1. Constructing the measure. Choosing one of the three vertices a,, i = 
{1,2,3} consider the unique ellipse Ei which: a) passes through ai and b) has 
dj , au as its foci. The constructed probability measure Mi is supported on the 
elliptic domain Ei bounded by Ei. We need the following notion. 

Given two distinct points ct\ ^ 02 on C define the arcsine measure wr ai)a2 i of 
the interval [01,02] as the measure supported on [ai,aa] and whose density at a 
point t G [ai,a 2 ] equals —7— — L_ _. 

To describe the measure Mi consider the family of straight lines parallel to the 
tangent line to the ellipse Ei at a^. Take now the family $i of intervals obtained 
by intersection of the latter straight lines with the elliptic domain Ei. Denote by 
— Vi the vector connecting a% with its opposite point on Ei, i.e. draw the straight 
line through a, and the center of Ei till it hits Ei again and take the difference 
of the latter and the former points. (One can easily check that if we introduce 
a new variable Zi = z — a, and express Q(z) = zf + Vizf + WiZi then the above 
vector will be exactly — v% in the expression for Q(z) which explains our notation.) 
Now parameterize the above family of the intervals by their middle points using 
the formula —Vi9 2 , 8 E [0, 1]. Consider the family fig of arcsine measures of these 
intervals. Finally the required measure Mi is obtained by the averaging of fig w.r.t. 
parameter 9, i.e. Mi = figd9, see Fig. [5p). 



4 



B. SHAPIRO AND M. TATER 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 -0.5 0.0 0.5 1.0 



Figure 2. a) The measure /z and the three ellipses Ei,E 2 ,E 3 for 
Q{z) = z(z — 1) (z + \ — i). b) The measure /i, ellipse E\, and 
several straight segments belonging to the family $1. 

Now we can finally formulate the main results of this paper. 

Theorem 2. // the measure fi in Conjecture [7] exists then each of the measures 
Mi, i € {1,2,3} have /i as its inverse balayage, i.e. fi and Mi have the same 
logarithmic potential (or the same Cauchy transform) outside the ellipse Ei and 
the support of /i is contained inside the support of Mi . 

By definition the Cauchy transform C v {£) and the logarithmic potential pot v (z) 
of a (complex-valued) measure v supported in C are given by: 

Cu{z)= ! ^ and pot v {z)= f log \z - £|di/(0- 

About the properties of the Cauchy trasform and the logarithmic potential of a 
measure consult e.g. j3]. 

Remark 1. Theorem [2] is so far a conditional statement. For technical reasons 
complete proofs of the existence, uniqueness and several other properties of fi are 
postponed until [12"] . 

Denote by CQ i (z) the Cauchy transform of the measure Mj, i = 1,2,3. The next 
result shows that each Cauchy transform Cq 1 (z) satisfies outside the elliptic domain 
Ei the following nice linear non-homogeneous second order differential equation 
(similar to the one obtained earlier in [3]). 

Theorem 3. The Cauchy transforms Cq^z) of the measures Mi, i = 1,2,3 de- 
fined in Theorem [H satisfy outside the ellipses Ei one and the same linear non- 
homogeneous differential equation: 

Q(z)C^(z) + Q'(z)C' Qt (z) + ^C Qi (z) + = o. (3) 
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2. Proof of Theorem [D 



Proof of Theorem [H It essentially follows from the stronger version of the main 
result of [7] which we present below. First we express the polynomial Sp n (X) as 
the characteristic polynomial of a certain matrix. In order to make this matrix 
tridiagonal we assume as above that the root a,; is placed at the origin. In order to 
simplify the notation we drop the index i assuming that z is already the appropriate 
coordinate. Set 



Q(z) = z 3 + vz 1 + 



Consider the operator 

T = (z 3 + 



d 2 

wz)—it + (az 2 
az z 



(3z- 



d 

l] Tz 



t (z-A), 



where v,w,a, /3, 7 are fixed coefficients of Q{z) and P{z) respectively and ni \ 
are variables. Assuming that S(z) = uqz 71 + u\z n ~ x + . . . + u n with undetermined 
coefficients m, < i < n, and in order to solve the Heine- Stieltjes problem described 
in the introduction we will be looking for the values of 9 n , A and u,;, < i < n, such 
that T(S(z)) = 0. Note that T(S(z)) is in general a polynomial of degree n + 1 
whose leading coefficient equals uq [n(n — 1) + an — 9 n ] . To get a non-trivial solution 
we therefore set 

On = n { n ~ 1 + a). 

Straightforward computations show that the coefficients of the successive powers 
z n , . . . , z° in T(S(z)) can be expressed in the form of a matrix product M n U, 
where U = (uo, Ui, 
matrix 



M n 



/A - £„,i 

In. 2 





. , u n ) T and M n is the following tridiagonal (n + 1) x (n + 1) 



a n ,2 

A — in,2 
7n,3 





A — £n,3 






«n.4 







A Cn,n 
7n,n+l 



with 



— i)(n — £ + 1) + (3(n — i + 1) 
(n — £)(n — i + 1) + a(n — i + 1) 



7n,i 



- 1, 



iu(n - i + l)(n - z + 2) + 7(71 - i + 2) 



IG{1,... 

ie{2,.. 
*e{2,. 





Ot n ,n+l 
A — Cn,n+1 



« + !}, 
.,»+!}. 



(4) 



A similar matrix can be found in [S] and also in [16j . The matrix M n depends 
linearly on the indeterminate A which appears only on its main diagonal. Obviously 
if the linear homogeneous system M n U = is to have a nontrivial solution U = 
(uq, u\, u n ) T the determinant of M n has to vanish. This gives the required 
polynomial equation 

S Pn {\) = det(M„) = 0. 

The sequence of polynomials {Sp n (X)} ne z + does not seem to satisfy any rea- 
sonable recurrence relation. In order to overcome this difficulty and to be able to 
use the technique of 3-term recurrence relations with variable coefficients (which 
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is applicable since M n is tridiagonal) we extend the above polynomial sequence by 
introducing an additional parameter. Namely, define 

Sp„,i(X) = det M„,i, ie{l,...,n + l}, 

where M n> i is the upper i x i principal submatrix of M n . One can easily check (see, 
e -g-j El P- 20]) that the following 3-term relation holds 

Sp„,i{\) = (A - €n,i)Sp n ,i-l(\) ~ 1pn,iSp n ,i-2(X), i G {1, . . . , n + 1}, (5) 

where is as in Q and 

i>n,i = ^n.iln.i, 1 6 {2, ...,«+ 1}. (6) 

Here we use the (standard) initial conditions Sp n ft(X) = 1, Sp n ,-i{\) = 0. It is 
well-known that if all £ n ,i's are real and all ip n /s are positive then the polynomials 
Sp n ,i(X), i G {0, ...,n+ 1}, form a finite sequence of orthogonal polynomials. 
In particular, all their roots are real. In our case however these coefficients are 
complex. To complete the proof of Theorem [5] we state the following generalization 
of Theorem 1.4] which translated in our notation claims the following. 

Theorem 4 (A. Kuijlaars - W. Van Assche). If there exist two continuous functions 
£(t) and iJj(t), t e [0, 1], such that 

lim &,n = f(r), lim tpi,„ = ^(r), We [0,1], 

i/(n-\-l)— >t z/(n+l)— >t 

i/ien the asymptotic root-counting measure p of the polynomial sequence {Sp n (A)} n ez_L 
{Spn,n+i{ty}n<£Z, + (if it exists) and the average M of the acsine measures given by 



; [C(r)-2 v /^M,C(r)+2y^)]' 

/iawe £/ie same logarithmic potential outside the union of their supports. 

Recall that for a pair of distinct complex number a± ^ 02 the arcsine measure 
W[ ai . a2 ] is the measure supported on [0:1,0:2] and whose density at a point t e 
[oi,02l equals — , 1 

Remark 2. Although Theorem [3] is not explicitly stated in [7] it is very similar and 
its proof is completely parallel to that of Theorem 1.4 from this paper. 

From the explicit formulas for and ip n ,i (see (0]) and ©) one easily gets 
£(r)= lim £ i>n = - v (l - r ) 2 , 

i/(n+l) — >t 

V(r) = lim ^,„ = - W (1-(1-t) 2 )(1-t) 2 . 

i/(«+1)-»t 

Notice that the above limits are independent of the coefficients a, (3, 7 of the poly- 
nomial P(z). □ 

Lemma 1. The parametric curve Y given in the above notation by the formula 
£(t) ± 2^ip(r), t e [0, 1] is the ellipse passing through the origin and given in 
coordinates x = Re(z), y = Im(z) by the equation 

aux 2 + 2a i2 xy + a 22 y 2 + 2ai 3 x + 2a 23 y = (7) 

where 

a n =C 2 +4D 2 , a 12 = —(AC + 4BD), a 22 =A 2 +AB 2 , 
a 13 = 2D(BC- AD), a 23 = —2B(BC - AD) 

and A = -Re(v), B = -Im(u),C = -Im(v), D = Re(u). 
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Proof. We express the functions £ and tp as 



-v(l — t) 2 = — vff 2 = —v sin 2 ip, 
-w(l - (1 - t) 2 )(1 - r 2 ) = -w(l 




• 2 2 

-tysin ipcos <p, 



(8) 



where re [0,1], 9:= 1-t6 [0,1], and sin^ := 0, p G [0,tt/2]. Then 

£(r) ± 2"Ji)j(t) = —v sin 2 (/? ± y/—w sin 2y>. 

Thus the curve T C C is given by the parametrization T(<p) = — v sin 2 tp ± 
y/—w sin 2<p, where v,w,z G C and G [0,tt/2]. Set to = u 2 , so that yj—w = hi. 
Then T has the form: 

T(ip) — — usin 2 tp ± iusin2c/3 = (— Re(v) — ilm(v)) sin 2 ±i(i?e(it) + i Jm(w)) sin2<^. 
We, therefore, get the following system for its real and imaginary parts: 

!x(p) = A sin 2 if + Bsin2ip 
y{<p) = C sin 2 tp + D sin 2tp. 

Here A = -Re(v),B = -Im(u),C = -Im(v),D = Re(u) and p G [— 7r/2,7r/2] 
since T is 7r-periodic. 

To show that T is an ellipse passing through the origin and satisfying |J2J) sub- 
stitute ((9]) into the expression anx 2 (tp) + 2ai2a;(</5)y(y) + a22y 2 ( ( fi) + 2ai3x(y)) + 
2a232/(v)i where the coefficients defined in the statement of Lemma [T] Sim- 

ple calculations then show then that the latter expression vanishes identically, i.e. 
for all values of ip. 

To prove that {7J describes a real ellipse (and not some other real affine quadric) 
consider the determinant 



A 



an 
an 
ai3 



012 

0-22 
C123 



Ol3 
023 





-4{BC - AD) 4 



It is well-known that if A is negative then we have a real ellipse (A > corresponds 
to an imaginary ellipse, i.e. an empty set of solutions). Thus unless BC — AD = 
(which describes the situation with all three roots of Q(z) being collinear) then T 
is a real ellipse. To find its semiaxes a and b we calculate the following quantities: 



6 := 



an O12 

012 a22 



4{BC~ ADf 



i ■= an + a 22 = A 2 + C 2 + 4(B 2 + C 2 ). 



It is known that the roots Ai j2 of the characteristic equation A 2 — iX + 6 = are 
equal to 2a and 26 2 , (in particular, both need to be positive) where a, b are the 
semiaxes of the ellipse under consideration. We arrive therefore at 



46 



46 



and Vi 2 - 46 = yJ{{A - 2D) 2 + (C + 2B) 2 )((A + 2D) 2 + (C - 2B) 2 ). For the sake 
of completeness the eccentricity c of our ellipse can be expressed as 



46 



\ y{(A - 2D) 2 + (C + 2B) 2 )((A + 2D) 2 + (C - 2B) 2 ). 



□ 



Lemma 2. The foci of the ellipse coincide with the two roots of the polynomial 
Q(z) different from the origin. 
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Proof. The coordinates of the centre c = (x c ,y c ) of our ellipse satisfy: 
aux c + a 12 y c + a i3 = fn A C 



a 12 x t + a 22 y c + 023 = OJ 2 2 

Recalling that Q(z) = z(z 2 + vz + w) ~ z(z 2 + vz + u 2 ) we need to show that the 
coordinates {xf,yj) of the foci / of T satisfy the equation: 



— v ± y/v 2 — Au 2 
/ = x f + iy f = - . 

To do this we express them through A, B, C, D. First, we see that Re(v 2 — Au 2 ) 
A 2 + AB 2 - C 2 - AD 2 . Using the relation: 



where r = ^/^ 2 + rj 2 we get 



and 



r = y/(Re(v 2 - Au 2 )) 2 + (Im(v 2 - 4u 2 )) 2 = 4c 2 = 
y/((A - 2D) 2 + (C + 2B) 2 )((A + 2D) 2 + (C - 2B) 2 ) 

Xf = ^± ^=^/4c 2 + (A 2 + AB 2 - C 2 - 4D 2 ) 
C . 1 



y f = — ± — y/Ac 2 - (A 2 + AB 2 - C 2 - AD 2 ) 

Straightforward calculation shows that the centre c and the foci /i and f 2 lie on 
the same line given by the equation: 

4c 2 -M 2 +4B 2 -C 2 -4L> 2 ) / A\ C 
y = 2(AC + ABD) L { X -2) + 2- 

Finally we check that the spacing between the centre and either focus equals 
to the eccentricity c which settles the lemma. This follows, for example, from the 
expression for the coordinates of the intersection points between (fit))) and the circle 
(x-x c ) 2 + (y-y c ) 2 = c 2 . □ 

3. Proof of Theorem 02 

We start with the following integral representation of the required Cauchy trasform. 

Lemma 3. The Cauchy transform Cq(z) of the measure Mq associated with the 
root of the polynomial Q{z) = z(z 2 + vz + w) at the origin is given by 

Cq(z) = f -= d ° (11) 

Jo yV - 4 ^) 614 + ( 2vz + ^w)6 2 + z 2 

Proof. Indeed, recall that the Cauchy transform C[ Ql Q , 2 ]of the arcsine measure 
^[ai,a 2 ] of the interval [0:1,0:2] equals 

c 1 

<-["i,«2] — / w r- 

yfz - ai)(z - a 2 ) 

The measure Mq is obtained by the averaging of the family of arcsine measures, 
namely 

M o = / wr t , . n cm.,.,, nnzftdr, 



l [C(r)-2^M^(T)+2 V / ^)]' 

where f (r) = -v(l - r) 2 = -v8 2 , ip{r) = -w(l - (1 - r) 2 )(l - r) 2 = -w{l - 6 2 )9 2 
and 9 = 1 — r. Since the Cauchy transform of the average of a family of measures 
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equals the average of the family of their Cauchy transforms one gets after obvious 
simplifications: 

C {z) 



o (z-(C(r)-2 v A^))(z-(C(r) + 2 v ^(7)) 
d6 

o yf{v 2 - 4w)6» 4 + (2vz + Aw)6 2 + z 2 ' 

□ 



3.1. Special case. We first provide the proof of Theorem [3] for a specific case 
Q(z) = z(4z 2 — 1) where the calculations are somewhat simpler and then address 
the general case. By Lemma [3] the Cauchy transform Cq(z) of the measure Mq 
associate with the root of Q(z) at the origin is then given by the integral 



C (z) 



dO 



Jo V6> 4 - e 2 + z 2 

We want to find a differential equation satisfied by Cq(z) w.r.t. the variable z. 
Unfortunately, we do not know how to do it directly and our proof requires a number 
of intricate variable changes and manipulations. We first change t = 29 2 — 1 and 
consider 

1 f 1 dt 



cw-jbM-^y^ ^ , (i2) 

where s := Az 2 — 1. Introduce now a family of functions I v (s) indexed by v > 
and defined by: 

T , . If 1 t v dt 
Iu{s) 



V2J-1 VtTTVt 2 + a ' 
Lemma 4. For v > the following three relations are satisfied: 



dL+2 l r dl v 
as 2 os 

^ + « = -i 7 ° + vrTi' (14 » 



Proof. Relation (flHj) can be proved directly: 
dh +2 1 f 1 t u+2 dt 



ds 2V2 J-i ^/tTl{t 2 + s) 3 / 2 

(t u+2 +st u )dt 1 f 1 st v dt 



2V2 J -l Vt+T(t 2 + Sf/ 2 2V2 J-! Vt+l(t 2 + S) 3/2 

l T dl v 
2 os 

Relation (jT4j) is easy to verify by integration by parts. Indeed, 

d_ 1 f 1 t + 1 tdt 1 1 

ds {2+ l} ~ 2V2j- 1 VtTT(t 2 + s)V 2 ~2VTT7s 4 °- 

Similarly, by integration by parts one gets: 
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a, , 1 r t + i tdt _3 I 

□ 

Now, we express 812/ ds from (fT5|) . substitute it in (|14[) . and single out dli/ds: 

^ = S 97 + I /0+ 27fTl- (16) 
Adding flU) with (JTBJ) and reducing dh/ds, dl 2 /ds with the help of (O we obtain: 



Through (fTOj) we get: 



4(l + s )^A = / + /i. 



.Wo 



Ji = 4s(l + s)— - + s/ + 2Vs+T. 
as 



Differentiating both sides of the latter relation w.r.t. s and using (|16|) again we 
obtain the required linear non-homogeneous differential equation satisfied by Iq(s): 

16,(1 + 16(1 + 2 S )^ + 3/o = - 7 =. (17) 

In order to recover the required equation ^ for Cq(z) we have to change s back 
to z. Using straightforward relations 

dl 1 dC , a 2 7 1 / <9 2 C aC 



and 



z- 



9s 8z 9z 9s 2 64z 3 \ dz 2 dz 
we obtain after some obvious simplifications the equation: 

z(4^ 2 - 1)^ + (12z 2 - 1)^ + 3zC (z) + 1 = 
dz 2 dz 

which coincides with ^ for Q(z) = z(4z 2 — 1). Thus our special case of Theorem[3] 
is settled. 

Notice also that (fT7|) can be solved explicitly. The general solution of the cor- 
responding linear homogeneous equation is an arbitrary linear combination of a 
complete elliptic integral of the first kind 2/1 (s) and of an associated Lcgendre func- 
tion of the second kind 2/2(5) given by: 




TTyT+S \ 2VT+i 

= Q_ 1/4 (l + 2s), 

here lK(x) and Q(x) are the complete elliptic integral and the associated Legcndrc 
function of the second kind respective. The general solution to (|T7|) depends on 
two arbitrary constants 61,62 and is given by: 

^ ^ f 9 ds f g ds 

lo = 612/2 + 622/2 + 2/2 / 2/1-TT77 ~ Vi / 2/2- 



f 2 W aL J a \f 2 W 7 
where = f 2 (s) = 16s(l + s), VK(s) = yi(«)^ a («) - 2/2 How- 
ever, we need its particular solution and thus have to determine the corresponding 
particular values of 61 , 62 . (To find them we evaluated the integral (fT2|) for two dif- 
ferent values of s. Moreover, analyzing the polynomial (t+l)(t 2 +s) = t a +t 2 + st+s, 
we observed that it is positive on [—1,1] for s > and that (fT2"|) is divergent when 
s = 0.) 
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The next figure compares the appropriate solution of (|17[) giving Iq(s) with the 
values of Jo(s) calculated numerically using the integral (jTSJ) for a number of values 
of s (which are shown by dots below). 




Figure 3. The graph of Io{s) obtained from (jTTJ) and its numerical 
values obtained from (fT2| . 



3.2. General case. The scheme of this proof is exactly the same as in the above 
special case but calculations are somewhat messier. Assuming that v 2 — Aw ^ 
we need to find a differential equation satisfied by the integral We change 

variables as follows: 



s = — 16w 



z + vz ■ 



v + 2z 
v 2 — Aw ' 



a = v — Aw 



(v 2 — Aw) 2 ' 

and denote Iq(s,u,o) = Cq(z). (Here as above we assume that v and w are some 
fixed complex numbers.) It also helps to change the variable 9 in (jTTJ) by using 
29 2 = t + 1 and then we finally get 



C (z) = J (s) 



1 



dt 



'2a J -i % /F+T v /(t + u) 2 + s 
As above we introduce a family of functions I v (s), s > given by the formula: 

us) : =^r ( '+ u)y * . 

v 7 ^ J -i VtTTVtt + u ) 2 + s 

Analogously to Lemma fj] one can prove the next statement. 
Lemma 5. The following relations are valid for J„(s), s > 0: 



dl v 



-2 



dl v 
as 



ds 



h) = u— -J 

as A 



d 



(h - h) = (u 2 



2u 



dh 



u-1 
-J-i) 



l) 2 + s 



(18) 
(19) 
(20) 



Now, we use (fT8"f for expressing 9J2 /9s and then we single out dl\ / ds from (p~9l) : 

(21) 



9s 



9J 



la 



1 



4(1 -u) 2^5(1- u) V("+l) 2 + s 
Adding Q19|) and (j2"U|) . employing (fT5|) again, and using (|210 we get the relation 



9Jn 



(« - l)/r = -4s(s + (u - lD^ - sJ - 



2(s-u 2 + l) 



I) 5 



(22) 
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Eventually, taking the derivative of the both sides of the latter equation w.r.t s 
and using (|21|) again we finally get a linear differential equation in the variable s 
satisfied by Iq(s): 



16 S (s+(u-lf) d ' 



^2^ dI 



2 s + (u + l)(5u+l) 

— -t- + 1o{2s+(u— 1) ) — — \~di -i — p= , , = — 

ds 2 ds V2 y/{ u + 1) 2 + s 



0. (23) 



In order to get an equation for Cq(z) w.r.t. the variable z, we use: 



and 



d 2 C _ d 2 s dip 
dz 2 dz 2 ds 



dCo 

dz 

ds 
Or. 



ds dlo 
dz ds 



du dip 
dz du 



d 2 I ds du d 2 I 



With the help of %h = 2I X and 



dCo 

dz 



ds 
dz 



du 

du 



dz 1 



dlo 
ds 



ds 2 dz dz dsdu 

we obtain 

du Iq du 



du 



d 2 I 
du 2 



dzA{l-u) dz ^(i_ u ) N /( M + i)2 + s 



Now, we get 

dlo 

ds 



(v 2 - Aw)(vz + 2w) dC v(v 2 - Aw) 



16wz 



dz 



32wz 



Co 



v(v 2 — Aw) 
32wz(v + z) 



(24) 



Further, we use 



d 2 I _ 4 gJp 
du 2 ds 



As 



d 2 I 
ds 2 



and 

d 2 Ip 
dsdu 



«(.-l>f£ + 



3s + A{u - l) 2 dl 



3/n 



3s + 5u 2 + 6u + 1 



2s(u-l) 9s 8s(u-l) 4VS(w-l)s(s + (w-l) 2 ) 3 / 2 ' 



We can now express d 2 Co/dz 2 through d 2 Io/ds 2 , dlo/ds, and Iq as follows: 



<9z 2 




4s 



3s + 4(u- l) 5 




s(u - 1) 
ds \ ( du 

dz) As(u-\) ' \d^J [dz) 4y/Z(u- l)s(s+ (u- l) 2 )3/2' 



This leads to: 



d 2 C 256wz 2 d 2 I 

dz 2 (v 2 — Aw) 3 ds 2 



16 



Aw 2 (w + z 2 ) + Avwz(w + 2z 2 ) + v 2 w(w + 5z 2 ) + v 3 (2wz — z 3 ) dl 



(v 2 — Aw) 2 (2w + vz)(w + z(v + z)) 
3v(v + 2z) 



A{2w + vz)(w + z(v + z)Y 
v(v + 2z){2,v A + 8v 3 z - 2Avwz - Aw(2w + 3z 2 ) + v 2 {bz 2 - 8w)) 
A(v 2 — Aw)(v + z) 3 (2w + vz)(w + z(v + z)) 
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From the latter equation and ([2"4"|) we finally get: 

d 2 I Q (v 2 - 4w) 3 d 2 C (v 2 - 4w) 2 d dC Q 

ds 2 256u>z 2 dz 2 256w 2 z 3 (w + z(v + z)) dz 

v(v 2 - Aw) 2 c 2 „ v(v 2 - 4w) 2 c 3 

~ 1024wz 2 (w + z{v + z)) ~ 1024w 2 z 3 (v + z) 3 (w + z(v + zj) ' 

where 

ci =4w 2 (w + z 2 ) + Avwz(w + 2z 2 ) + v 2 w(w + 5z 2 ) + v 3 (2wz - z 3 ), 
c 2 =8vwz + v 2 (w - 2z 2 ) + Aw(w + 4z 2 )), 

c 3 =w 4 (w - 2z 2 ) + 12vwz(w + 3z 2 ) + 4wz 2 (3w + 4z 2 ) + v 3 (8wz - Az 3 )+ 
+ v 2 (4w 2 + 27wz 2 -2z 4 ). 

Plugging these formulae into (|23p we arrive at: 

4z(z 2 + vz + w)Cq(z) + 4(3z 2 + 2vz + w)C (z) + (3z + v)C (z) + 1 = 0, 

which can be equivalcntly expressed as 

Q(z)Co(z) + Q'(z)C (z) + Q"(z)C (z)/8 + Q'"(z)/24 = 0, 

with Q(z) = 4:z(z 2 + vz + w). (Notice that the multiplication of Q(z) by a non- 
vanishing constant is irrelevant in our considerations.) □ 

4. Final remarks 

It is very tempting to extend the methods and results of the present paper to 
the case of the 'generalized' Hcun equations which are of the form 

{ Qk+l{z) Jz^ + Qk{z) d^ + - + Q2{z) Tz + V{Z) ) S{Z) = °' 

where degQfc+i(z) = k + 1 and dcgQ;(z) < i for i = 2, 3, k. As in the introduc- 
tion for each positive (and sufficiently large) integer n there exist n + 1 polynomials 
V(z) counted with appropriate multiplicities such that for each of these V(z) the 
above equation has a polynomial solution S(z) of degree n. Thus one can define 
the corresponding spectral polynomials and study the asymptotics of their root- 
counting measures. Large scale numerical experiments support the following. 

Conjecture 3. For any 'generalized' Heun equation the sequence {fJ. n } of the root- 
counting measures of its spectral polynomials converges to a probability measure \x 
supported on a curvilinear planar tree located inside ConvQ k+1 and whose leaves 
(i.e. vertices of valency 1) is the set of all roots of Qk+i(z), see Fig. [7} Moreover, 
the limiting measure /i depends only on Qk+i(z), i.e. is independent of the other 
coefficient of the equation. 

We finish our paper with the following problem. 

Problem 1. Under the assumption that the latter conjecture holds (which is very 
likely) is it true that the Cauchy transform C M of the limiting root-counting measure 
fi satisfies a linear ode of the form: 

Qk+iWP (z)+a x Q' k+1 {z)C^- 1 \z)+a x Ql +1 (z)C^ (z)+...+a k+1 Q ( j ^ 1 1) (z) = 0, 

where ai,...,ak are some universal constants, i.e. independent of Qk+i(z) (but 
maybe dependent on the order k of the operator). 
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Figure 4. The measure fj, for the operator Q(z)4^ with Q(z) = 
(z - 5)(z - I)(z - 4 + I)(z - 2 + 31) (2 - 3 - 27). * 
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